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We present a determination of the pion-nucleon ( nN ) cr-term a n N based on the Cheng-Dashen 
low-energy theorem (LET), taking advantage of the recent high-precision data from pionic atoms to 
pin down the nN scattering lengths as well as of constraints from analyticity, unitarity, and crossing 
symmetry in the form of Roy-Steiner equations to perform the extrapolation to the Cheng-Dashen 
point in a reliable manner. With isospin-violating corrections included both in the scattering lengths 
and the LET, we obtain ovjv = (59.1 ± 1.9 ± 3.0) MeV = (59.1 ± 3.5) MeV, where the first error 
refers to uncertainties in the nN amplitude and the second to the LET. Consequences for the scalar 
nucleon couplings relevant for the direct detection of dark matter are discussed. 
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INTRODUCTION 


The nN cr-term measures the amount of the nucleon 
mass that is generated by the two lightest quarks. Since 
the dominant contribution originates from the energy 
content of the gluon field, due to the trace anomaly of the 
QCD energy-momentum tensor, the nucleon mass would 
only change moderately if the quark masses were turned 
off. Thus, UkN encodes information on the explicit break¬ 
ing of chiral symmetry and constitutes one of the funda¬ 
mental low-energy parameters of QCD. In recent years, a 
precise determination of the cr-term has become increas¬ 
ingly urgent, given its relation to the scalar couplings of 
the nucleon that are prerequisite for a consistent inter¬ 
pretation of direct-detection dark matter searches [lj-H] • 

Traditionally, information on avjv has been in¬ 
ferred from nN scattering by means of the Cheng- 
Dashen LET Hi that relates the Born-term-subtracted 
isoscalar amplitude D + at the Cheng-Dashen point s = 
u = m 2 N , t = 2M 2 , to the scalar form factor of the nu¬ 
cleon cr(t) evaluated at t = 2M 2 (precise definitions be¬ 
low) . The application of the LET thus requires two main 
ingredients: the analytic continuation of the isoscalar nN 
amplitude into the unphysical region, and the correction 
due to the finite momentum transfer in <r(2M 2 ). The first 
task has been addressed by extrapolating partial-wave 
analyses (PWAs) from the physical region to the Cheng- 
Dashen point by means of dispersion relations @-0], hi 
particular, in [3,LLQJ a formalism was developed to express 
the result of the extrapolation in terms of threshold pa¬ 
rameters for nN scattering. Similarly, the scalar form 


factor requires a dispersive reconstruction to account for 
the strong 7T7 t rescatteringin the isospin-0 S -wave Ill- 
Based on the PWA from Q, a value ct^n ~ 45 MeV 
was inferred in [10]. This result was later challenged 
by a new PWA [121 ] . leading to a much larger value of 
VnN = (64 ± 8) MeV, although based on the same for¬ 
malism. In fact, the discrepancy could be traced back, 
to about equal parts, to different input for the isoscalar 
nN scattering length, the nN coupling constant, and nN 
partial waves for the evaluation of the dispersive inte¬ 
grals. 

A second strategy that has been pursued relies on Chi¬ 
ral Perturbation Theory (ChPT) to perform the extrap¬ 
olation to the Cheng-Dashen point. However, to deter¬ 
mine low-energy constants still input for the nN phase 
shifts is required, so that the outcome of the ChPT anal¬ 
yses tends to support the value of a^N corresponding to 
the PWA used as input EMU. Moreover, it has been 
questioned whether the chiral representation is at all ac¬ 
curate enough to permit a reliable extrapolation to the 
Cheng- Dashen point 15]. For a detailed comparison to 
results obtained in lattice QCD, we refer to |16| . 

In this Letter, we combine two new sources of informa¬ 
tion on nN scattering that have become available over 
the last years. First, the measurement of level shifts and 
decay widths in pionic atoms 17-fl9l] has led to a_preci- 
sion determination of the nN scattering lengths EflEH. 
Second, a system of Roy-Steiner (RS) equations has been 
developed [22| that combines general constraints on the 
nN scattering amplitude imposed by analyticity, unitar¬ 
ity, and crossing symmetry. The construction proceeds 
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similarly to Roy equations for 7T7 t scattering 1231 ]. where 
the solution for the low-energy phase shifts can be param¬ 
eterized in terms of the 5-wave scattering lengths 24]. 
In the case of nN scattering, the construction and solu¬ 
tion is complicated by the presence of the crossed channel 


7T7T —» NN, cf. [25l - l27| . as well as the increased number of 


relevant partial waves. While partial results have already 
been presented in [H, Ei - dio! ], here we use the complete 
solution of the RS system to obtain, in combination with 
the scattering-length constraints from pionic atoms, a 
precision determination of the nN cr-term. In particular, 
at this level of accuracy the impact of isospin-violating 
(IV) corrections cannot be ignored, as demonstrated by 
the isoscalar nN scattering length 20. EE EH EH , so that 
revisiting the Cheng Dashen LET becomes mandatory. 


where A# subsumes higher-order corrections in the chi¬ 
ral expansion. Corrections to the LET have been investi¬ 
gated systematically in SU (2) ChPT, with the result that 
An is very small: non-analytic terms are absent at full 
one-loop order BEl , so that the dominant corrections 
are expected to scale as ~ 1 MeV. Indeed, 

estimating the low-energy constants based on resonance 
exchange, one obtains [33| 

|A fl |<2MeV, (8) 

an estimate which we will adopt in the following. 

In practice, 0 is usually rewritten as 

ovat = T, d + A fl — A a — Ar, (9) 


where 


CHENG-DASHEN LOW-ENERGY THEOREM 

We start by stating the precise formulation of the 
Cheng- Dashen LET [2, 5j. In the isospin limit, the scat¬ 
tering amplitude for the process 

n a ( q ) + N(p)^n\q') + N(p / ), (1) 

with pion isospin labels a, b and Mandelstam variables 

s = (p + q) 2 , t={p'-p) 2 , u = (p — q') 2 , (2) 


A ct = a{2M 2 ) - a nN , A D = D+( 0,2 M 2 ) - E d , 

^ d = F 2 (4 0 + 2M 2 d+ 1 ). (10) 

Here, F\ = 92.2 MeV [3] denotes the pion decay con¬ 
stant and the subthreshold coefficients are defined via 
the expansion 

OO 

o + {v,t)= y d+y m t n . (ii) 

n,m—0 


fulfilling s + t + u = 2m 2 N + 2 M 2 , can be expressed as 

T ba {v,t) = 8 ba T + (v, t) + \[r b ,T a ]T-{v,t), 

T\v,t) = u(p')|-D 7 (i/,t) - j^^£ 7 (z/,t) ju(p), (3) 

where v = (s — u)/(4toat), I = ± refers to 

isoscalar/isovector amplitudes, uin and M^ to the nu¬ 
cleon and pion mass, r a denotes isospin Pauli matrices, 
and we normalize spinors as uu = 1. Amplitudes A 13 
with definite s-channel isospin I s are 

(£)-(!-.)(£)■ • 4efD>B) ' (4) 


The LET involves the Born-term-subtracted amplitude 

.2/1 1 


D + (v,t ) = D + (v,t) — 


9 

m N 


vg 


L N 


V N 


(5) 


where g is the nN coupling constant. The scalar form 
factor of the nucleon is defined as the matrix element 


a(t) = {N (p')\rh(uu + dd)\N(p)), m =— ( 

( 6 ) 

with up- and down-quark masses m u and m d , momentum 
transfer t = (p' —p) 2 , and <r(0) = cr^N- The LET then 
states that 


D+(0,2M 2 )=a(2M 2 ) + A R , (7) 


Although individually sizable due to strong nn rescat¬ 
tering, the difference Ac — A CT was shown to be small 
in |ll]. Here, we use the updated value |35, 


A D ~A a = (-1.8 ± 0.2) MeV, 


( 12 ) 


which incorporates modern input for 7T7 t phase shifts, ef¬ 
fects from KK intermediate states, and the uncertainties 
due to nN parameters. 

As alluded to above, the isoscalar channel is known 
to be sensitive to IV corrections. For this reason, we 
now derive a version of the LET that takes the dominant 
IV effects into account. First, we define the cr-term as 
the average value of proton and neutron scalar-current 
matrix elements {N S {p, n}) 


&7 tN - 


(Jr) I (Jn 


&n = (N\m(uu +dd)\N), (13) 


where, up to third order in the chiral expansion, one finds 

(Jp — CTr) 


37]. Next, we identify the isoscalar amplitudes 
everywhere with the average of the n±p —> n ± p charge 
channels 




X p =n(X. 


rr+ T> ^-ir r 


2 V' ~ 7T'p — yTV'p I VV 7T p—>TC 


-p), (14) 


for X € {D, doo, cioi,...}. The motivation for doing so 
is two-fold: first, the n±p charge channels dominate the 
nN data base, so that this scenario is closest to the one 
considered in PWAs. Second, the uncertainties in the nN 
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scattering lengths are smallest if one works in the phys¬ 
ical, not the isospin basis [20j, 21 1. As a consequence, 
we identify the nucleon and pion mass with the masses 
of the proton and the charged pion, respectively. We 
also assume that the radiative corrections applied in the 
PWAs remove the dominant effects, and therefore we con¬ 
sider all quantities to be virtual-photon subtracted. In 
this scenario, the leading IV corrections are generated 
by the mass difference between charged and neutral pion 
Att = M 2 — M 2 0 . For the scalar form factor one finds [37] 

=a p (2Ml)-a p (15) 

= + 9 128 nF* ( ~~ 7 + 2 ^ l0S (1 + ’ 

where gA denotes the axial charge of the nucleon. Simi¬ 
larly, the IV corrections to Ap can be extracted from |38j] 

A p D = F 2 {d p ( 0,2 Ml) - d p 00 - 2MX} (16) 

= + 9\M„A / 4^1 (1 + V 2 )). 

384t 1F 2 T 256t tF 2 \ SV ’) 

Taking everything together, we obtain 

o'ttN = F 2 (cIq 0 ± 2M 2 <ig 1 ) + A-d — Ac- — A# 

Slg^M^A^ e 2 2 / \ 

+ 2567 tF 2 +Y K ^ h + f2 > 

= F 2 (d p 0 + 2MXi) + (1-2 ± 3.0) MeV. (17) 

In 1171) we also included the leading corrections due to 
virtual photons, encoded in the low-energy constants /1 
and f‘ 2 . The latter can be determined from the proton- 
neutron mass difference [ 13 ], fi = (—0.97±0.38) GeV -1 


Since the errors are dominated by different sources—IV 

1/2 
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for the former we use the estimate |/i| < 1.4 GeV -1 
Hot ]. The single largest correction is generated by A v , an 
upward shift of 3.4 MeV. Such large IV corrections have 
already been observed in the case of the 7 tN scattering 
lengths 31. ±2 . 


PIONIC ATOMS 

Pionic hydrogen (71 -H) and deuterium (7 tD), a n~ and 
a proton/deuteron bound by electromagnetism, provide 
access to nN physics due to the imprint of strong inter¬ 
actions in the energy spectrum. The shift of the ground- 
state energy level in nH and ttD, as well as the width of 
the 7 riJ ground state, probe three different combinations 
of 7 tN scattering lengths. The input quantities relevant 
for the RS equations are the s-channel-isospin scattering 
lengths ciq!,, defined in terms of the n±p charge channels. 
Updating the analysis of Hip to account for the new 
value of the 7 riJ level shift [19j and subtracting virtual- 
photon effects as detailed in J41J, we obtain 


a ,= (169.8 ± 2.0) x 10 ~ 3 M~\ 
a,= (-86.3 ± 1.8) x 10 


corrections in the case of a 0 (, and uncertainty in the ex¬ 
traction of the isoscalar combination for a^J 2 —the errors 
can be considered uncorrelated to a very good approxi¬ 
mation. Apart from their role in the solution of the RS 
equations, the scattering lengths are also a crucial ingre¬ 
dient in the determination of the irN coupling constant 
via the Goldberger-Miyazawa-Oehme sum rule [42]. In¬ 
deed, if the scattering lengths from :IJ are used, one 
recovers the value g 2 /( 47r) = 14.3, whereas CHI leads to 
5 2 /(47 t) = 13.7±0.2 [ 20 , P- This result, to be adopted in 
the following, stands in good agreement with more recent 
determinations from NN [13] and itN 441 scattering. 


ROY—STEINER EQUATIONS 

Roy equations [ 2 ^] for 7T7 t scattering, or RS equa¬ 
tions 2f| 25l- 271 for non-totally-crossing-symmetric pro¬ 
cesses, incorporate the constraints from analyticity, uni- 
tarity, and crossing symmetry in the form of dispersion 
relations for the partial waves. They can be shown to be 
rigorously valid in a certain kinematic region, in the case 
of irN scattering the upper limit is s m = (1.38 GeV) 2 jp . 
The integral contributions above s m as well as partial 
waves with l > Z m , with Z m the maximal angular momen¬ 
tum explicitly included in the calculation, are collected in 
the so-called driving terms, which need to be estimated 
from existing PWAs, as do inelastic contributions below 
s m . The free parameters of the approach are subtraction 
constants, which, in the case of 7T7 t scattering, can be di¬ 
rectly identified with the scattering lengths 2 4j, while for 
the solution of the nN system it is more convenient to 
relate them to subthreshold parameters instead. The re¬ 
sulting system of coupled integral equations corresponds 
to a self-consistency condition for the low-energy phase 
shifts, whose mathematical properties were investigated 
in detail in 145;]. Following j24j, we pursue the follow¬ 
ing solution strategy: the phase shifts are parameterized 
in a convenient way with a few parameters each, which 
are matched to input partial waves above s m in a smooth 
way. To measure the degree to which the RS are fulfilled, 
a \ 2 -likc function is defined according to 


x 2 = ^ fi^-FlfkW 

i,i.,±j= 1 v 


R emWj) 


(19) 


(18) 


where {Wj} denotes a set of points between threshold 
and ^/s m , are the s-channel partial waves with isospin 
I s , orbital angular momentum Z, and total angular mo¬ 
mentum j = l ± 1/2 = Z±, and F[/y[] the right-hand 
side of the RS equations. We take l m = 1, N = 25 (dis¬ 
tributed equidistantly), and choose the number of sub¬ 
traction constants in such a way as to match the num¬ 
ber of degrees of freedom predicted by the mathematical 
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FIG. 1: Phase shifts 5j £ of the s-channel partial waves in de¬ 
grees, obtained from the solution of the RS equations. The 
dashed line indicates our central solution, the bands the un¬ 
certainty estimate. The partial waves are labeled by the spec¬ 
troscopic notation I/ 2 / s 2 j- 


properties of the Roy equations [45]. It should be stressed 
that the form of the RS equations only reduces to that of 
Roy equations once the f-channel is solved, see 22]. In 
the solution of the RS equations we minimize HU with 
respect to the subtraction constants (identified with sub¬ 
threshold parameters) and the parameters describing the 
low-energy phase shifts, while imposing (fT51) as additional 
constraints. 


earizing around the central values m , we find 
E d = (57.9 ± 0.9) MeV + ^ cj.Aa^, 

Is 

Cl /2 = 0.24 MeV, c 3/2 = 0.89 MeV, (20) 

where Aagf, measures the deviation from (1181) in units 
of 10Already in this linearized form, one re¬ 
covers T,d from 0 if the scattering lengths from [5, [S] 
are inserted, while the modern input produces E^ = 
(57.9 ± 1.9) MeV (this also indicates that the 5-wave 
phase shifts from f| need to be amended close to 
threshold). Moreover, the difference to [p2j can be traced 
back to the P- wave scattering volume a^, , which needs 
to be known extremely accurately due to its large weight 
in the formalism of [0, 0. Once the RS equations are 
solved, the threshold parameters can be calculated from 
sum rules, and indeed we find that the result for af + is 
slightly lower than the value used in [l3|, which already 
suffices to explain the difference 0- The main impact 
of the RS equations in the u-term determination thus 
amounts to eliminating the need for independent input 
for a ^_, . In total, our result for the er-term becomes 

VtvN = (59.1 ± 3.5) MeV. (21) 

Although already 4.2 MeV are due to new corrections to 
the LET (thereof 3.0 MeV from isospin breaking), we do 
observe a significant increase compared to 0. As illus¬ 
trated by m, this effect can be immediately traced back 
to our modern knowledge of the ttN scattering lengths as 
extracted from pionic atoms. By combining this informa¬ 
tion with the constraints from RS equations, the er-term 
can be determined to a remarkable accuracy. 


We performed a number of checks as regards the sensi¬ 
tivity of our solution to the input quantities: the number 
of grid points N as well as the number of parameters 
used in the description of the partial waves were varied, 
the s- and f-channel partial waves in the driving terms 
truncated at different Z max = 4,5 and J max = 2,3, the 
matching conditions at s m as well as the s-channel partial 
waves evaluated from different PWAs, and the sensitivity 
to the precise definition of the y 2 -function was investi¬ 
gated. In addition, to stabilize the fit we imposed sum 
rules for the higher subthreshold parameters. The solu¬ 
tion for the s-channel partial waves, expressed in terms 
of the phase shifts and including uncertainty estimates 
from these systematic studies as well as the uncertain¬ 
ties in the scattering lengths and the coupling constant, 
is shown in Fig. [T| A more detailed account of our RS 
solution will be given in 0. 

Apart from low-energy phase shifts, the RS solution 
provides a consistent set of subthreshold parameters. In 
particular, this allows us to pin down E^ in accord with 
both the RS and the scattering-length constraints. Lin- 


SCALAR NUCLEON COUPLINGS 

The existence of a weakly-interacting massive particle 
(WIMP), one of the most promising dark-matter can¬ 
didates, could be established in direct-detection experi¬ 
ments, which are sensitive to the recoil of the WIMP scat¬ 
tering off nuclei (see 0] for a review). The interpreta¬ 
tion of these searches relies on the couplings of the WIMP 
to nucleons, according to its quantum numbers. A pre¬ 
cise determination of ct^n therefore has immediate con¬ 
sequences for the scalar channel, since, as it was shown 
in [3|, the scalar couplings of the nucleon to q = u, d, 

m N fq = ( N\m q qq\N), (22) 

follow once cr^N is determined, with all further correc¬ 
tions taken into account within SU( 2) ChPT. Taking 
m u /md = 0.46 ± 0.03 from [47[, we obtain 

fP = (20.8 ± 1.5) x 10" 3 , fP = (41.1 ± 2.8) x 10“ 3 , 

/" = (18.9 ± 1.4) x 10" 3 , fd = (45.1 ± 2.7) x 10” 3 . 

(23) 






















5 


In addition, we quote our result for 

E = l+\{fu+fd+fs) = 0.305±0.009, (24) 

q=u,...,t 

avera ged over proton and neutron, and with /j v taken 
from |48| (in principle, the strangeness coupling follows 
from the cr-term by means of SU( 3) considerations, but 
the uncertainties are too large to compete with recent 
lattice determinations). This particular combination of 
scalar coefficients becomes relevant in the context of 
Higgs-mediated interactions, not only in direct detection, 
but also in Higgs-induced lepton flavor violation [49| ■ In 
particular in (1231) the uncertainties have been appreciably 
reduced, thanks to the precise knowledge of cr^N inferred 
from our RS equation analysis of nN scattering. 
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